Superradiant light scattering and grating formation in cold atomic vapours 

R. Bonifacio", B.W.J. McNeil 6 , N. Piovella a , G.R.M. Robb b 
a Dipartimento di Fisica, Universita Degli Studi di Milano, 
INFN and INFM, Sezione di Milano, Via Celoria 16, Milano 1-20133, Italy 
b Department of Physics and Applied Physics, University of Strathclyde, 
Glasgow G4 ONG, Scotland, UK 

A semi-classical theory of coherent light scattering from an elongated sample of cold atoms exposed 
to an off-resonant laser beam is presented. The model, which is a direct extension of that of the 
collective atomic recoil laser (CARL), describes the emission of two superradiant pulses along the 
sample's major axis simultaneous with the formation of a bidimensional atomic grating inside the 
sample. It provides a simple physical picture of the recent observation of collective light scattering 
from a Bose-Einstein condensate [S. Inouye et al., Science 285, 571 (1999)]. In addition, the model 
provides an analytical description of the temporal evolution of the scattered light intensity which 
shows good quantitative agreement with the experimental results of Inouye et. al. 

PACS numbers: 42.50.Fx, 42.50.Vk, 03.75.Fi 



I. INTRODUCTION 

Recent experiments by S. Inouye et al. at MIT |IJ] have demonstrated the formation of atomic matter waves in a 
cigar-shaped Bose-Einstein Condensate (BEC) pumped by an off-resonant laser beam, together with highly directional 
scattering of light along the major axis of the condensate. This emission has been interpreted as superradiant Rayleigh 
scattering, and some theoretical work describing this experiment has been recently published [^]||. In particular, the 
work of Moore and Meystre [|| describes the Rayleigh scattering in a BEC using a model which extends the Collective 
Atomic Recoil Laser (CARL) model originally proposed by Bonifacio et al. j|-[7j to include a quantum-mechanical 
description of the centre-of-mass motion of the atoms in the condensate |8]|§]. The conclusions of ref. || were that 
the original CARL theory, which treats the atomic centre-of-mass motion classically, fails when the temperature of 
the atomic sample is below the recoil temperature Tr — Huj r /kB, where u> r — h\q\ 2 /2m is the recoil frequency, m is 
the atomic mass, q = k — k s is the difference between the pump and scattered wavevectors and is Boltzmann's 
constant. However, the cubic dispersion relation derived in ref. || reduces to that of the original semiclassical CARL 
model for large atomic densities: More specifically, the quantum corrections to the classical motion are negligibly 
small when the CARL parameter p in the free electron laser (FEL) limit roughly interpreted as the average 

number of photons scattered per atom, is greater than one. This suggests that a fully quantum-mechanical description 
of the atomic centre-of-mass motion may not be necessary in order to describe the main experimental results of ref. 
0, i.e. the temporal evolution of the scattered light intensity and the spatial grating in the condensate. We are 
aware that a semiclassical theory is necessarily limited in its description of the radiation statistics and the quantum 
degenerate nature of the condensate, which require a full quantum analysis. Nevertheless, we consider the semiclassical 
approach useful in order to give an intuitive description of the physical mechanism underlying the observed effects. 
We stress however that in spite of its simplicity, the semiclassical model produces good quantitative agreement with 
the experimental results of ref. Q. 

II. MODEL 

The model described is bidimensional and semiclassical. We represent the cigar-shaped atomic sample as an ellipsoid 
with length L and diameter W, where L ^> W as shown in Fig|l]. The sample is exposed to a classical plane wave 
radiation electric field Eo(y,t) — x£oe ik ^ v ~ c:t ' > + c.c, polarised along the x axis and incident along the axis y, with 
£ o real and constant and where k = lu/c We assume that the scattered radiation consists of two radiation pulses 
propagating along the z axis, with electric fields polarised as the incident field: 

E(z, t) = x[£x(z, ty k{z - ct ^ + £ 2 (z, t)e~ tk{z+ct ^ + c.c], (1) 

where £^2(2, t) are slowly varying complex amplitudes. The dominance of the scattering along the ±z axis over that 
in other directions is due to the geometry of the atomic sample. 
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The atomic sample is described as a collisionless gas of atoms, each with two internal energy levels. The internal 
evolution of each atom is described by the density matrix elements p mn (to, n = 1,2) for the lower, (1), and upper, 
(2), levels. The off-diagonal elements pi 2 = p\\ describe the dipole moment induced by the radiation fields via the 
relation d — xp{p\ 2 + c.c), where p, is the dipole matrix element. The diagonal elements p\\ and P22 describe the 
probability of an atom being in the lower or in the upper level, respectively. The off-diagonal elements may be 
described conveniently as a sum of three polarisation waves: 



P12 = S e lk{y - Ct) + S ie ik( - Z ~ ct) + S 2 e- lk{z+ct) . 



(2) 



The dipole moment of each atom contributes to the macroscopic polarisation of the atomic sample described by 
P = n(x)d, where n(x) is the atomic density. This polarisation is a source for the radiation field via Maxwell's wave 
equation which yields, in the usual Slowly Varying Envelope Approximation (SVEA), 
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where we have neglected the terms proportional to e ± 2luJt . We assume that the atomic sample can be described as a 
collection of N point particles with positions Xj, so that n(x) =£f=i<5 (3) ( x ~ Xj). Multiplying both sides by e Tlkz 
and integrating over the z axis from z — Az/2 to z + Az/2, where Az — A/2, Eq.(||) yields 
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where the upper sign corresponds to the first subscript and we have assumed the field amplitudes £^2 are spatially 
slowly varying over Az. Assuming also that £1.2 are independent of x and y, we can integrate on the plane (x, y) over 
the section A = nW 2 /4 of the condensate, so that Eq.(||) becomes 
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where n = N/AAz is the average density and (..) = (1/N) YljLi(--)j- 

In this model the atomic centre-of mass motion is treated classically, with each atom described as a point particle 
with a given position and momentum. The radiation fields drive the centre-of-mass motion of the atoms via the force 
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Neglecting the fast-varying temporal terms, the equations for the atomic velocity components are: 
dv„. 
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We assume that the detuning S = lo — uj a between the optical fields and the atomic resonance is much larger than the 
natural linewidth of the atomic transition, 7, so that the atoms always remain in their lower internal energy states 
(pu ~ 1 and p 22 m 0). Moreover, assuming that the scattering time scale is much longer than the relaxation time 7 , 
we can adiabatically eliminate the atomic polarisations, i.e. Sk = *(/V^)£fc/(7 + ~ ^fc/2#, where Qk = 2fj£j~/h, 
|f2fe| is the Rabi frequency for the field k and k = 0, 1, 2. With these approximations, Eqs.(|g) and (0) yield: 
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It is seen that the interference between the pump and the scattered fields forms two bidimensional periodic potentials 
2(2/1 z) cx |£o£i,2| cos[k(z =F y) i 0i,2] m the plane (y, z), where <fii^ 2 are the phases of the complex amplitudes £ x 2 . 
A weaker ID potential V^(z) cx |£i£2| cos[2fcz + <f>x — (f> 2 ] forms along the z axis due to the interference of the two 
counterpropagating scattered fields. If the pump intensity is large enough, we can assume £q 3> £1,2 and neglect the 
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ponderomotive potential V3. Then, Eqs.(||),(^) and (^) can be conveniently written in the following dimensionless 
form @]: 
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where 8i 2 = k(z =F y), pi^ — (m/hkp)(v z =p t^) and = —2i(e /hunp) 1 / 2 £i2 are scaled atomic position, atomic 
momentum and field amplitude variables respectively. The dimensionless time and space coordinates, t — cu r pt and 
z = ui r pz/c, are scaled in terms of the collective recoil bandwidth, pui r , where u> r — hk 2 /m is the single-atom recoil 
frequency and p = (Qo/2d) 2 / 3 (ujp 2 n/eouJ 2 h) 1 / 3 is the dimensionless CARL parameter At t = 0, the atoms are 

assumed to be randomly distributed in position and have zero momentum, and the amplitudes of the scattered fields 
are set to zero. 



III. ANALYSIS 



In this simple model the two scattered fields are uncoupled and symmetric. For each field (1,2) individually, Eqsinn)- 
( |l2| ) are formally identical to those which describe pulse propagation in a high gain free electron laser (FEL) fjJTlt 
is already known that they admit a self-similar solution of the form Ai^{z,t) = ±zA(u), where u = \/pi[(I =p z) and 
A(u) is the solution of a set of ordinary differential equations |LG(| . This self-similar solution describes the superradiant 
emission of radiation pulses whose duration decreases in proportion to the fourth root of the peak intensity. The pulse 
shape can be approximated by a hyperbolic secant function, followed by some non-linear 'ringing', similar to that 
which occurs in superfluorescence from inverted two-level atoms pd| . 

A simpler model can be obtained by approximating the spatial derivative in the field equation ( [l2|) by a damping 
term [H i.e. 
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where k = c/2ui r pL and L/c is the transit time of the photon along the major axis of the condensate. In this 
approximation, the finite interaction time due to the escape of radiation from the atomic sample is represented by 
an incoherent decay of the field amplitude in the sample at a rate c/2L, half the inverse of the radiation 'lifetime' in 
the atomic sample. A more general treatment where the radiation is scattered in a direction making an angle ip with 
respect to the z axis should give k ki (c/2uj r p)[\ simp\/W + | cos ip\/L] p^ |. As L ^> W, the radiation is least strongly 
damped along the major axis of the sample. 

An approximate solution to eqs.(|l0|),(il"|) and (|l|) can be found assuming k 3> 1 and adiabatically eliminating 
the field variables, i.e. A X 2 ~ K _1 (exp[=Fi6' 1 2 ])- In this limit, the rate of change of the average scaled momentum 
is (d/dt)(pi t 2) — z f2k|Ai.2| 2 - A third-order analysis of the equations in the mean-field limit (i.e. with radiation 
propagation modelled by the damping term) gives the following approximate solution: 
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where tjy = — %/2/cln (|6o|/v2) is the delay time of the peak and 60 = (exp[=Fi6'i j 2(i = 0)]) is the initial bunching, 
which can be assumed to be ~ \j\pN for a condensate of N atoms. In the linear regime the exponential gain is 



G = Wrpy/2/K = (3j/S) ^ (2I N/miu){\ 2 /A), 
whereas the peak value of the scattered intensity is 

I P eak = {j/6) 2 [(3/2ir)(\ 2 /A)N} 2 I , 



(16) 
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where Iq = 2ceo|£o| is the pump intensity, A is the cross sectional area of the condensate and 7 = p, k /Girheo is the 
natural decay rate of the atomic transition. 
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IV. COMPARISON WITH THE MIT EXPERIMENT 



In the MIT experiment, a sodium BEC was exposed to a single off-resonant laser pulse red-detuned by S/2ir = 1.7 
GHz from the 35*1/2 ~ * 3/3/2 transition, with A = 0.589/im and natural width 7 = 0.31 x 10 8 s _1 . The recoil frequency 
is uj r = 3 x 10 5 s" 1 . We assume that the condensate had a diameter of 20 p m and a length of 200/im, approximately 
N = 5 x 10 5 atoms participate in the emission of a scattered radiation pulse. The dimensionless parameters are 
p 44 x I^ 3 and k « 5.5 x 10 4 x I Q 1 ^ 3 , where Iq is the pump intensity in mW/cm 2 . As Jo > 1 and consequently 
p > 1, the results of ref. |s) indicate that quantum effects due to atomic diffraction should be negligibly small for 
this experiment, even though T <C Tr. As k > 1, the atoms emit two superradiant pulses along the major axis of 
the condensate. The gain is approximately G w 82 x \/Io, where G is given in ms , whereas the peak occurs after 
a time to — ln(2A)/G w (170/\/7o)/xs, in good agreement with the measured values of ref. m. Furthermore, from 
Eq.(|l5j) the modulus of the average atomic velocity is v = {Tik/m)p\{pi t 2)\ ~ (X/2ir)G s» 0.7\/iocrn s _1 . Fig . [2] shows 
the temporal evolution of the main peak of the scattered intensity, as given by the approximate formula (|1J), for the 
parameters of the MIT experiment and three different values of the incident intensity, 3.8 mW/cm 2 (solid line), 2.4 
mW/cm 2 (dashed line) and 1.4 mW/cm 2 (dotted line). 

In addition to the temporal evolution of the scattered radiation pulses, there are other predictions of this semiclas- 
sical model which are consistent with the results of the MIT experiment: Firstly, superradiance is observable only 
if the Doppler broadening of the atomic resonance is sufficiently small that atsR <C 1, where a is the rms spread of 
the gaussian spectral distribution and tsR ~ 1/G is the superradiant time The observed spectral width of the 
Bragg resonance of the BEC of approximately 5 kHz Q] (corresponding to a velocity spread of few mm s~l) yields 

ctsR ~ 0.16 x I Q X ^ 2 . We observe that, using a = fc(^T /m) 1 / 2 , a temperature of only IpK (approximately the BEC 
transition temperature for the MIT experiment) would increase the frequency spread by a factor of 15, enough to 
destroy the superradiant emission. This explains why superradiant emission was observed only at the extremely low 
temperatures below the threshold for Bose- Einstein condensation Q. Secondly, superradiant emission parallel and 
antiparallel to the z axis induces an average atomic velocity v\ t 2 ~ (G/k)[y^fz\, respectively at 45° degree with respect 
to the negative (positive) direction of the z axis, as observed in the MIT experiment. We assume the existence of two 
distinct families of atoms interacting with the two independent superradiant pulses A\ and Ai- However different 
orders of atomic velocity, i.e. v m ,n — mv\ + nv^, with m,n integers, have also been observed. More precisely, the 
orders (2,0), (—1,1), (0,2), (2,1) and (1,2) other than the usual (1,0) and (0,1), have been clearly observed in the 
experiment after increasing the exposure time to the laser source and letting the atomic cloud expand ballistically 
after the interaction ^] . The formation of this momentum grating can be explained as a sequential superradiant scat- 
tering process in which the atoms emit m pulses along the positive z, and n pulses along the negative z axis, acquiring 
a total recoil velocity v m>n . The extremely narrow resonance line allows the atoms to emit up to three sequential 
superradiantly scattered pulses before at ~ 1, which is consistent with the observation of the atomic momentum 
distribution. 



V. CONCLUSIONS 

In conclusion, we have presented a semiclassical model describing the superradiant Rayleigh scattering from a Bose- 
Einstein condensate observed in Ref. jlj] . The model is much simpler than those previously used to explain the results 
of Q| as the atomic centre-of-mass motion is treated classically. The evolution of the scattered intensity and the atomic 
motion due to recoil as calculated from this simple model are in quantitative agreement with the experimental results. 
The fact that quantum centre-of-mass effects such as atomic diffraction are negligible is a consequence of the high 
density of the condensate. In our model the BEC is essentially described as a collisionless Doppler-free atomic gas. 
The results presented here suggest that together with its high density, the most important property of the condensate 
with regard to superradiant light scattering is its very low temperature rather than its quantum degenerate nature. 
In this respect the situation is similar to that of ultraslow propagation of light in a BEC [ p"5| . Subsequent observations 
of ultraslow propagation in a hot vapour |l(| demonstrated that the significant property of the BEC was that it was 
a Doppler-free optical medium rather than a quantum degenerate one. 
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FIG. 1. The geometry of the scattering esperiment. The filled elipsoid, representing the atomic condensate with dimensions 
L and W , is illuminated with a single off-resonant laser beam of electric field Eq polarised along the x axis and propagating 
along the ±z axis. The geometry favours the emission of the oppositely directed superradiant pulses E\ and E2 along the 
major axis of the condensate. 

FIG. 2. Temporal evolution of the main peak of the scattered intensity as given by the approximate formula (14), for the 
parameters of the MIT experiment and three different values of the incident intensity, 3.8 (solid line), 2.4 (dashed line) and 1.4 
(dotted line) mW/cm 2 . 
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